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Abstract 

Working within the framework of the covariant perturbation theory, we obtain the 
coincidence hmit of the heat kernel of an elliptic second order differential operator 
that is applicable to a large class of quantum field theories. The basis of tensor 
invariants of the curvatures of a gravity and gauge field background, to the second 
order, is derived, and the form factors acting on them are obtained in two integral 
representations. The results are verified by the functional trace operation, by the 
short proper time (Schwinger-DeWitt) expansions, as well as by the computation of 
the Green function for the two-dimensional scalar field model. 



1 Introduction 



The heat kernel is a unique mathematical subject that represents a fundamental 
block for building any quantum field theory [1]. This subject was extensively studied 
in mathematics |i2j and in physics [3] . For the general review of the current state of 
research we recommend [1] . The present paper is a technical but important step in the 
development of one of the directions in heat kernel research, the covariant perturbation 
theory (CPT) proposed by A. Barvinsky and G. Vilkovisky [5]. In this paper we 
use notations and techniques of the preceding papers on the covariant perturbation 
theory [SI 13 E], however, to keep this work self-contained we include their basic 
formulas. This may be needed because the bulk of [H [9] is not published in journals 
or electronically. 

The setup and algorithms of CPT were outlined in [6] , whose main ideas are briefly 
repeated in this Introduction and at the beginning of the section O We consider 
quantum field theories in space-times with the metric tensors of constant positive 
signature, i.e., Euclidean space-timesll] The space-time dimension 2a; is arbitrary 
throughout the paper, except for the section [51 If S[(^] is a classical action of a field 
theory, then the induced differential operator is, 

^■-(V)*(^.«) = ^^SH. (1) 

The operator -F(V) acts on small fluctuations of an arbitrary set of fields (^'^{x). 
Throughout the paper, like in any other work on the covariant perturbation theory, 
we employ DeWitt's condensed notations [10]. In this convention the field index A 
stands for any set of discrete indices of tensor-spinor fields which allows us to describe 
a wide range of field models from scalar to non-Abelian gauge fields. That is it, all 
results presented here are valid for arbitrary background non-Abelian fields. Also, in 
DeWitt's notations the integrals over the spacetime are omitted except for the overall 
integral. 

The heat kernel K{s\x,y) is assumed to be a solution of the heat equation in the 
form, 

^^k{s\x,y) = F{y^)k{s\x,y), (2) 

where s is the proper time [T2l ITOt [TT] along the geodesic connecting two space-time 
points X and ?/, and operator -F(V) is acting on fields at point x. The hat notation 
in ([2]) and elsewhere indicates matrix valued quantities. The matrix conventions are 
1 = 5^Bi P = P^Bi etc. The matrix trace over the index set A is denoted by 
tr, e.g., trP = P^a- We restrict ourselves to quantum field models that can be 
characterized by the elliptic second order differential operator of the type called in 
physics literature a minimal second order operator p!3]. 



F(v) = ni + p-im, (3) 

o 



^The metric g^'^{x) is characterized by its Riemann and Ricci curvatures R'^.^^p 
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where the Laplace-Behrami operator (or Laplacian), 

□ = g^'V^V,, (4) 

is constructed in terms of the covariant derivative V^, which is specified by its com- 
mutator curvature 

[v^, v.]<^^ = (v^v. - v.v^)<^^ = n^B^.u^''. (5) 

The exphcit presence of the Ricci scalar term in the operator ([3]) rather than in the 
potential is related to conformal models studies in four dimensions [6] . The set of the 
field strengths 

7t^,, P (6) 

characterizing the background will be called the curvatures denoted by the sym- 
bol 3?. Manifolds under consideration are asymptotically fiat with trivial topology, 
specifically, their gauge and gravitational curvatures and potential ([6]) vanish at the 
space-time infinity [6]. The work [13J contains a detailed review of the classes of 
models which are associated with the operator ([3]) and, therefore, can be dealt with 
by the present formalism. 

The formal operator definition of the heat kernel is 

k{s\x, y) = exp [sF(V^)] 6{x, y). (7) 

Knowledge of the heat kernel with separated points K{s\x^y) would allow one to 
construct the covariant diagrammatic expansion to all loop orders [TOj [H] . Here we 
focus only on the coincidence limit of the heat kernel K{s\x, x). As a result, quantum 
field theory applications are limited to the one-loop approximation. In particular, the 
one-loop effective action is defined as [13], 

W = ^Tr InF - J S'^x5^^'^\x, x){. . .), (8) 

and it represents the generating function for the one-particle irreducible diagrams 
[15] . This fact is especially important for quantum gravity where the work with 
ordinary Feynman diagrams would be formidable. The second term in ([8]) stands for 
the contribution of the local functional measure [16]. As shown in [17] this contribu- 
tion always cancels the volume divergences. For the massless operators ([3]) the result 
of such cancellation is a subtraction of the zero curvature term in the heat kernel 
expansion. Then, the one-loop effective action W is given by the trace of the heat 
kernel 

-(Trir(s)-Trir(.)U=o), (9) 

Z JO s 

where Tr denotes the functional trace 

Trir(s) = j (f'^xiik{s\x,x). (10) 



W 
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From this point, we do not indicate the space-time dimension at the integral measure, 
i.e., dx = (P^ X. The knowledge of the heat kernel, in contrast to its trace ffTOj) . would 
provide one with the Green function or propagator G{x, y) defined by 

F(V^)G(x,y) = -i5(x,y), (11) 

due to the Schwinger equation, 

G{x,y)= / dsK{s\x,y). (12) 
Jo 

The ultraviolet divergences appear in quantum field theories as divergences of the 
proper time integrals (Q, f|T2l) at the lower limit. They are to be removed by the 
renormalization procedure [TOl [13] . For massive field theories, there is a mass factor 
g-sm ^Yie heat kernel (trace), which makes the integral converge at the upper limit 
s ^ oo. Early techniques (Schwinger-DeWitt expansion) could only reproduce the 
asymptotic behaviour of the heat kernel K{s) at s ^ 0. As a result, these proper 
time integrals for massless field theories would diverge at the upper limit. These 
infrared divergences due to the method of computation were irrelevant of a quantum 
field theory under study. This is clearly seen from the fact that, upon the proper time 
integration, the short time expansion of the heat kernel corresponds to the inverse 
large mass expansion of the effective action or the Green functions ^U\. Therefore, a 
method that would allow one to compute the heat kernel in the whole range of the 
proper time values was needed. Such a method, proposed by G.A. Vilkovisky [18] 
and called the covariant perturbation theory, was systematically developed in a series 
of works O |6[ [71 [9]. CPT is the covariant expansion of the heat kernel by orders 
of the curvatures. The CPT heat kernel contains an infinite number of derivatives, 
expressed as form factors, acting on the background field curvatures. Thus, it is a 
nonlocal expression, and the heat kernel form factors turn into the Green functions 
acting on the curvatures. For the heat kernel trace up to the third order in the 
curvatures it has the general form [191 [20] (a different method leading to the nonlocal 
effective action was developed in [22]): 

TtK{s) = fdxg'/'{x)tTll + s^ + s'j2Ms,a,,n,)^,^^{t) 

29 ^ 

+ X: F.(s, Di, Ds, □3)3fJi3fJ23fJ3(0 + 0[3?^] , (13) 
i=i ) 

where g is the determinant of the metric tensor, and fi and Fi are for analytic functions 
(form factors) of the dimensionless quantities 

e = -sn, (14) 

which act on tensor invariants constructed of the curvatures 3?. The index of an 
operator □ indicates the local curvature it is acting on, i.e., □2-R1-P2 = R{x){pP{x)). 
The expressions resulting from these operations, in the curly brackets in f[T3l) . are 
taken at the same spacetime integration point x. 
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The CPT calculations are carried out with accuracy 0[3fJ'^], meaning the derived 
expressions contain terms of up to n — 1 order in the curvatures explicitly. In con- 
trast to the short proper time (Schwinger-DeWitt) expansion of the heat kernel, the 
expression f[T^ is valid for any value of the proper time s. In general, as is shown in 
[B], the large time behavior of TyK{s) is 

TrK{s) oc s'^"^, s^oo, 3? ^ 0, (15) 

for all curvature orders, except the zeroth. This property guarantees the convergence 
of field theory integrals at the upper limit in space-time dimensions 2uj > 2. 

In this work we derive the coincidence limit of the heat kernel ([7]) up to the 
second order in curvatures. Even though the final result is the same, this task is 
accomplished by two different methods. One method is a direct application of the 
CPT formulas [6]. The other is the generating function method. The method of 
generating expressions for the Green functions first has been proposed for analysis of 
local divergences of the coincidence limits of the Green functions [Ml [23]. It reflects 
the fundamental feature of the effective action as the generating function of the one 
particle irreducible Green functions. In [21] it was shown that this method can be used 
to treat nonlocal curvature expansions as well. The work with heat kernels, instead of 
the Green functions, has the advantage that one is not restricted to specific space-time 
dimensions, and no divergences are present until the proper time integrals are done. 
Since the differential operator -F(V) depends on three independent background fields, 
the metric, the gauge field, and the potential, there are three possible variational 
equations [HI [23l |2l]. The one we are interested in is with respect to the potential 
P, 

k(s\x,x) = --I—TTK{s). (16) 
s 5P{x) 

We start from the trace of the heat kernel f|T3l) . The only important feature of the 
form factors in f|T3|) required here is that they are functions of the operator f|T^ and 
do not depend explicitly on the curvatures 3?. This is the feature of the variational 
method which always eliminates one curvature reducing an accuracy 0[3?"] by one 
order. In this case, the variation of form factors is not required. By this approach 
we can only obtain the terms valid to the second order in the curvatures of the form 



K[s\x,x) = — 



I i=l i=l 



(17) 

In the next section, we also derive this expression directly by the CPT algorithms 
and confirm the equivalence of two methods. 

An important element of the present study is use of computer symbolic manipu- 
lations. There are two different kinds of symbolic manipulations performed. One is 
purely algebraic manipulations of form factors which was performed with general pur- 
pose mathematics software Maple. The other is manipulations with tensors in order to 
analyze tensor invariants, but especially to compare these results with the Schwinger- 
DeWitt coefficients from other works. These manipulations were performed with help 
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of MathTensor ^26j and Ricci [27] packages, both of which work under Mathematica 
software. 

This paper represents a natural part of the series of works on the covariant per- 
turbation theory El El El E] , and it consists of six sections. Sections [2] to [3] contain 
main results, namely, the heat kernel coincidence limit up to the second order in the 
curvatures. We start from the basic perturbations expansion and proceed to the co- 
variant expressions. Its form factors are obtained in two integral representations. In 
section m we study the short proper time asymptotic behaviour of K{s) and compare 
these asymptotics with the known Schwinger-DeWitt coefficients; the large (late) time 
asymptotic is reproduced as well. Section [5] is different from others, since it deals with 
a particular class of field models, the Weyl invariant model in two dimensions. We 
compute the Green function for this model and show that the result coincides with 
the known closed form solution. In Summary we address possible applications of the 
obtained results for QFT models. 



2 Covariant perturbation theory for the heat ker- 
nel 

Direct computations according to the covariant perturbation theory rely on whose 
basic equations we reproduce here to make the present paper self-contained. We also 
emphasize that the work with form factors of the second order of the heat kernel is very 
much similar to the work with the third order in the heat kernel trace. Therefore, in 
this and the following sections we use notations and techniques developed in [3 El [9] . 

To set up the perturbation theory, splitting of the metric and the covariant deriva- 
tive into auxiliary parts and perturbations is introduced: 

gf^u ^ ^M^' + Zi/^'^, (18) 

Vf,^ = V^y^ + t^^. (19) 

The auxiliary metric and derivative are taken to be fiat, i.e., 

R^^uapig) = 0^ (20) 

[V^,VJ^ = 7^^,VP = 0. (21) 

Now there are three independent perturbations, one each for the metric, the con- 
nection and the potential 

f^, P'Ir^- (22) 
In terms of these perturbations the differential operator ([3]) is divided as follows, 

F(V) = 5 + V(V), (23) 

with the fiat Laplacian 

□ = r^^.^u (24) 
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and the perturbation 

Viy) = /i/^-V^V. + 2f^V^ + P - ^M, (25) 

where 

^ (^M^ + /^A'^)f ^ (26) 



(otherwise, raising and lowering indices are done with help of the auxiliary metric 
)• 

The heat kernel is expanded then in powers of the perturbation [121 E]: 

oo 

K{s) = Y: Knis) (27) 



n=0 



where Kn{s) is a term of the n-th power in the perturbations (122!) . Perturbative 
solution to ([2]) can be found iteratively and has the general form, 



rS rtrt ft2 — — — — 

Kn{s) = dtn dtn-i ... dhe'^'-'"^''Ve^'--'-'^°V . . . e'-''-''^°Ve''° , (28) 
Jo Jo Jo 

which is in fact merely the Dyson series [ISIEH]- For the simple case of fiat space-time 
without background gauge fields, the perturbation is just the potential term P. In 
this closed form of Kn{s) for any n can be easily written down [30] . 

An exact solution for the zeroth order of the heat kernel is known [TOl. 



Kois\x, y) = ^-l^^V4(^)^i/4(^),-^~^(^^ yy (29) 

Here Sq is the parallel transport operator along the geodesic connecting y to x, g is 
the determinant of the auxiliary metric g^'^, and a is the auxiliary world function (in 
the Cartesian coordinates it is simply (x — ?/)^/2). When calculated by the algorithms 
of [6], term Kn{s) is obtained in the generic form. 



1 



K4s\x, x) = -£—fl\x) I d-a 5(1 - Y: 



X exp 
With the notation 



1 



sfi„+i(ai,...a„+i|V*)J ^sC„(ai,...a„+i|xi) ^. (30) 



f ci"a(5(l-f:aO/(ai,---an|a;0 _ =(/)„, (31) 

Jai>0 Xi—x 



this can be rewritten 



1 " 

= TIT^V^'^' ^ ^'(e^^"-C7i)„. (32) 
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Here Qn{<^i, ■ ■ ■ ttn|V*) is an operator of the second order in V*, and the derivatives V* 
act on the perturbation with index i contained in C^. Each term in C*^(ai, . . . an\xi), 
where z = 1, . . . n, is a product of n perturbations fl22|) at the points xi, . . . x„ respec- 
tively, and the index i on a perturbation means that the perturbation is at the point 
Xi, e.g, Pi = -P(xi). After the action of V* is done, all points Xi are made coincident 
with the point x in (130!) . For the consistency with [3 [8], in the n-th order, we use 
the notation Qn+i instead of Qn- 

Let us now display the forms for Q and C, for n = 0, 



AO 
•-^0 



0, 



(33) 
(34) 



The results for n = 1 are 

f^2(ai,a2|V^) 
C'°(ai,a2|xi) 

Cl{ai,a2\xi) 

and for n = 2 



-hi. 



-R]l 
6 



P 



(35) 
(36) 

(37) 



(]3(ai, 02,03! V*) 
C2{ai,a2,a3\xi) 
C2{ai,a2,a3\xi) 



a2a3^i + aia302 + 0102^^ 



-hih2 H — 
4 2 

-Pih2 - 

~g,,{Dl + Dl) + 2~g,^{Dl - D^))f?/ir 



^R^h2-\{g^,DlDl + ~g^^DlD^ 



+ 4~g,^DlDl)K''h^ 



(A 



-Rii) (P2 - -i?2i 

6 ^ 6 



-2{Dl + Dl){P,-^R,i)Ti^ 

-2{DlDlDl + DlDlDl)rihi^'^ 

+ ADlDlT^f^2 + DlDlDMht^hfl 



(38) 
(39) 



(40) 



(41) 



here 



h^'-g 



(42) 
(43) 
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Dl = aiVj-aaVi, 



(44) 



Dl = a^Vl-a.Vl- (45) 

The variational method was apphed directly to the results for the trace of the heat 
kernel in terms of perturbations derived in [8j, and it resulted in the same expression 
above. 

Finally, one should make the series ( l30l) manifestly covariant. This means we need 
to replace the perturbations and the auxiliary metric and derivative (fT8|) - (fT9|) by the 
respective covariant curvatures ([6]), metric and derivative ([5]). A covariant expansion 
to the second order in the curvatures involves the first order in the perturbations 
as well. The curvature expansions of the perturbations (122]) can be obtained from 
eqs. fl20l) - fl21l) . Their general solutions are integral equations which were solved by 
iterations in terms of R^u and IZ^^ in [6j, Eqs. (4.28), (4.29). Alternatively, one can 
use the generating function method f|T6|) applied to Tri^(s) of f|T3|) . It is easy to 
observe that only thirteen of the 29 tensor structures of the trace of the heat kernel 
presented in [T^] contain matrix P and thus contribute to the heat kernel itself. 

Whether we obtain the covariant heat kernel by the generating function method 
or by direct computations, the result is 

k{s\x,x) = j^g^'^{l + s(g^{-sU)P + g^{-sU)Rl) 

5 

W E ^.(-sni, -sDa, -s□3)^^l^^2H 

i=l 
11 

i=6 

+ s^Gi2{-sU^, -sDs, -sn3)3?i3?2[12] + 0[^^]}. (46) 
Here 3?i3?2[*] with i = 1 to 12 are quadratic structures: 



5Ri3?2[l] = PiP2, 


(47) 




(48) 


5Ri3fJ2[3] = A^22, 


(49) 


3fJi3?2[4] = /2i/?2i, 


(50) 


gfJigfj2[5] = i?fi?2M4, 


(51) 


3f^l3f^2[6] = v^7^^v"7^2a., 


(52) 


3f^l3f^2[7] = [v,A,v^7^f], 


(53) 
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3fJi3fJ2[8] = V^V.Ai?r, (54) 
5Ri5J2[9] = V„i?i^,V^i?ri, (55) 
3?i3?2[10] = V^V,i?ii?ri, (56) 

3fJi3fJ2[ll] = VaVpRi.uV^V'Rfi. (57) 

There is an additional quadratic structure linear in TZ^u 

3?i3?2[12] = V.7^^V^i^2, (58) 

which is separated from the others because it is absent in the final answer. Without 
gravity, the basis flTn) - flS51) reduces to only four non-vanishing curvature structures. 

3 Form factors of the nonlocal heat kernel 

The form factors gi[—sn) and Gi{—sOi^ —3^2, —sO^) are obtained as the integrals 
over parameters (l3T|) and, in this form, represented by two nonlocal kernels (cf., ( l35ll 
and (I38D): 

exp{saia20) (59) 

and 

exp(sfi), = 020^01 + aias02 + 0102^3. (60) 
The function ( !59l) appears in the combinations 

Am = e'"^'"^°'", (61) 
Bm = '—^ , (62) 

■J '—'to 

Umn = _□ ^ > (63) 

VTOn = -), (64) 

where (m, n = 1, 2, 3; m ^ n) and the indices of the form factors refer to the indices 
of the laplacians appearing in them. The function ( I6OI) appears in the combinations, 

e^^, e^^ - 1. (65) 

The coefficients of these functions are polynomials in a 's, Laplacians, and inverse 
Laplacians (Green functions). It is assumed that Di acts on the curvature 3fJi, ^2 
acts on the curvature 3?2, and 

□3 = (Vi + V2)'. (66) 
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The operators of form factors commute with the derivatives V's acting on the 
curvatures in fl56|) - fl57|) . because contributions from such commutations belong to 
the higher order 0[5R'^]. 

The first order form factors are formed of the ones in fll3l) as, 

9i{0 = 2/4(0, (67) 

92(0 = /3(0, (68) 

here and below ^ is the dimensionless operator- valued argument f|T^ . The second 
order form factors ^1(6,6,6), ^^i i = 1 to 11, are expressed via form factors 
Fj(^i,,^2,6) of trace of the heat kernel ffTSl) . tabulated in [S], in the following 

way. 



Gi(6,6,6) 


= ^1(6, 6, 6) + i^i(6, 6, 6) + 


6,6), 


(69) 


^2(6, 6, 6) 


= -^3(6, 6, 6), 




(70) 


^3(6, 6, 6) 


= i^6(6,6,e2) + i^6(e3,6,e2)-^(e3- 


-6)i^i7(6,6,6), 


(71) 


^4(6, 6, 6) 


= ^^4(6,6,6) + ^(6 -6 -6)^^15(6 


, 6, 6), 


(72) 


^5(6, 6, 6) 


= -^5(6, 6, 6), 




(73) 


G^6(6, 6, 6) 


= -^14(6, 6, 6), 




(74) 


^7(6, 6, 6) 


= —-^13(6,6,6), 




(75) 


^8(6, 6, 6) 


= i^i7(6,ei,e3) + i^i7(6,e3,ei). 




(76) 


^9(6, 6, 6) 


= -^16(6,6,6), 




(77) 


^10(6, 6, 6) 


= —-^15(6,6,6), 




(78) 


^11(6, 6, 6) 


= -^26(6,6,6), 




(79) 


^12(6, 6, 6) 


= -^32(6, 6, 6), 




(80) 



where some interchanges and cyclic substitutions of indices of the arguments 6 are 
made. 

It should be emphasized that these rules are not sensitive to the representation 
of form factors. When polynomials of a -parameters are present in (13T|) . we call it 
a -representation. In a -representation the first order form factors are 



91 = (A),, (81) 

* = I) A -b)^. (82) 

And the second order form factors admit the form 

Gi = (e^^^)3, (83) 

G2 = (2aia2e^^^)3, (84) 



^ \2n2 S 2 02 S 2 ^2 /2 
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- 2- 



1 e 



□2 S 



+ 



□ 



I I I lo 

+ — (2aiQ;2 - "2) + 7:-(«2 - 2aiQ;2) 



□ 



+ 
+ 



1 , 2 1x^2 1 i32 

□ iDs 2 ' Di ^, 



s I2 



1 e"^3 _ 1 



□ 102 



+ 



□3 / . 
\?)a\0L2 — 2q;2 ) 



1/2 

+ — (—- + 4a2 + 2ai — 2^2 — 6q;iQ;2 



J2 3 



+ 



(51 2^2 r, 
1 — «! H — q;2 H — q;iQ;2 — 0:2'' 
36 2 3 3 

— 'oaxa'i + 2q;iq;2^ + 3q;i^q;2^) 

□2 ^3 6 

- a\a2 + 2q!iq;2^) + 7=^f - ^Q;2 + q;2^ - \oi.xOL2 

□2^6 3 

+ 4q;iQ;2 — ^ol\ ol^ — laxa-f + 2q;i q;2 ) 



□.2 



□102 



(ai a2 ~ Q!iQ;2^ 



1 



l'-'2 



+ 



+ 



3 □, \ B: 



□ 



1 ' 2x^ °3 4 2; 

-ai ) + 7-7— (-«! ) 



□in2 2 



□1 2nin2/ s 



1 C*^3 

-2 h 



1 1 

— (2q;iQ;2) + — (2q;iQ;2) 



□ 102 
1 



-2Q;ia2) 



11-12 



(2Q;2)e' 



8^3 



+ 2 



1 U 
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1 



(4a2')e' 



2\ 8^3 



□2 S / 2 

1 



1 



8Q;ia2)- 



□2 s / 2 
1 



□1^2 

1 ^3 



2xA 



□1^2 



(-2«i )-f 



+2 



□ 102 S2/2 



□1^2 



(4q;2 - 12q;2')^ + 



1^1 
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The a-polynomial representation for the form factors fl83|) - fl94|) is not unique, i.e., 
contributions of some form factors Gi, that in fact identically vanish, can be present 
in the heat kernel. Due to this fact there is an extra quadratic structure ( 158|) which 
will be absent in the final answer for K{s). Thus, we need to proceed to the explicit 
representation introduced in p] [9]. This has been done again in two ways, the first 
is the treatment of (l8Tl) - (!9^ with the a -polynomial reduction procedure derived in 
[HI E] , see Appendix A. The second is the use of expressions (l67j) -( |79l) and the data 
for the form factors in the third curvature order of the heat kernel trace [8], [9]. Both 
methods result in the expression fH6|) for the heat kernel where the contribution of 
the structure fl58|l vanishes. The explicit form of the first order form factors is as 
follows 

giiO = fiO, (95) 

1 1 f _ 1 

^2(6) = 1^/(0 + 2^^^^- (96) 

The structure and complexity of the next order form factors are similar to those of the 
third order form factors in the heat kernel trace. The full table of these form factors 
can be found in Appendix B. These tables, in form suitable for computer symbolic 
manipulations, can be obtained online from Nuclear Physics B or from the author. 

The obvious check of the calculations above is the functional trace operation flTUl) . 
We can use the a -representation of the form factors for this consistency check. All 
total derivative terms should be discarded, and all second order form factors reduced 
to first order ones by identities like 

tr Jdxg^/\x)F{-sn^, -sB^, -sn.^)^^^2 = -tTjdx ^^/^(x)/(-sn2)3?i3?2 + 0[^% 

^ . (97) 

While first order terms collapse to a single local term P, the tables of the second 

order form factors reduce to five form factors of the trace of the heat kernel found in 

Therefore, while performing the functional trace operation ffTOl) after the generat- 
ing function variation (fT6j) . we start with the higher order of the heat kernel trace and 
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end up with the lower order one. This shows that there exists a connection between 
two consecutive orders in the curvature of the heat kernel trace. Namely, each lower 
order in the curvatures is completely defined by the higher order. That is it, 

/C„_i = tr y dxj^ JCn, (98) 

where 

/C„ = tr 1 dxg'/^ ^^(Vi, • • • , V„)3f?i • • • 3fJn(^) (99) 

i 

is a given order in the curvatures. 



4 Short and large proper time behaviour of the 
heat kernel 

Historically, a great deal of work on heat kernels has been devoted to the short time 
expansion also known as the Schwinger-DeWitt (or HaMiDeW) series p^Of [T3]: 

1 oo 

K{s\x, x) = -^—^g'/\x) J2 ^"an(x, x), S^O. (100) 

The Schwinger-DeWitt coefficients ) have been calculated explicitly for n = 

to 3 in many works [TOl |13[ [31] . Several results for n = 4 were obtained for special 
cases [3D1 E2|, as well as for a general set of fields [23 [S3]- [S3] gives even the 
fifth coefficient, but its expressions are written in the index free form, therefore, it 
would require a special study to transform it into a representation suitable for our 
comparative analysis. All d„(x,x) coefficients are local functions of the curvatures 
([6]). The short time expansion serves as a limiting case and as a consistency check 
for the results presented in the previous section. 

To find the short time asymptotic expansion (11001) . we need to know short time 
expansions of the form factors in K{s\x,x). Generally, this procedure is very much 
similar to the one performed for the trace of the heat kernel in [H [20] . It is a simple 
exercise to expand the basic form factors f[59l) - f[65l) . By inserting these expansions 
into the table of form factors, we obtain the Schwinger-DeWitt series. Since a„(x,a;) 
are local functions of the background fields, the nonlocal denominators 1/ A cancel in 
the short time expansions of form factors, as expected. Still, the tree nonlocal terms 
l/D are present due to the Riemann tensor expansion in terms of the Ricci tensor 
[5]. We omit here these nonlocal expressions for the Schwinger-DeWitt coefficients. 

Technically, the tree nonlocal expressions appear due to the Bianchi identity [6] 

R'"^^'^ = ^(yV^R'"^ - VV'R''^ - V^V'^i?''" + V'V'^i?^") + O[i?2].(101) 

Since we are working in the second order in the curvatures, we also need the nonlocal 
expansion for the square of the Riemann tensor: 

RfiuafBR'^"'^^ = 'iRfifals'^ a — Rul3 + 0[R^] 
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1 



1 1 



li-'2 



4nin2 
+ 0[i?3]. 



□2" + Da" + 2nin2 - 2nin3 - 2n2n3)i?r^r 



(102) 



The true local representation of the derived short time expansion for the nonlocal 
heat kernel can be restored, if we find an irreducible basis of quadratic tensor invari- 
ants containing the Riemann tensor, for each of the coefficients. It is easy to do in 
the second order, as the basis for a2(x, x) consists of just one structure. 



(103) 



Since there are several derivatives eliminating the Riemann tensor with help of the 
Bianchi identity in the third and fourth coefficients, the basis for a3(x, x) may contain 
only 

n{R^,^pR^'^-P)^ R^^^f.V^V^R"^, (104) 

and for aJx, x), 



u\R^,^pR>^-'-P)^ □(i?^,,^V^V"i?'^^), □(/?^,,^nV'^V"i?^^). 



(105) 



We construct now the local a„(x, x) with help of the invariants fll03p - fll05p and other 
acceptable tensor invariants fl47p - (|37|) . mixed with a proper number of the Laplacians 
□ . These expressions are substituted into constructed local ) with unknown 

numerical coefficients. Equating the resulting quantities to the short time expan- 
sion in its nonlocal representation we arrive at a unique solution for the numerical 
coefficients to be found. 
The final result is. 



a2[x,x) 



a3(x,x) 



6 



+ 
1 

60 

+ 



—uRi + -pp + J-7^„^7^^^ 

180 2 12 



180 
1 



Rii,,R^'^ 

180 ^ 



-U^Rl 



1260 

U{PP) + —upp + —PUP 
24 ^ ^ 24 24 

90 ^ ^ ' 180 



360 ^ ' 



180 

— UPR + —PUR 
360 360 



-TZ^'^un 



ioU 



1 r, 
60 



13a 



90 



+ 



5040 



U{R>"'R^ 



1260 



1890 
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15 



^ -RnR+ —^n{RR) 



7560 15120 

(x,x) = —n^P + ^—n^Rl 
^ ' ^ 840 15120 



+ 0[^% (107) 



+ —n(pnp) + —npnp + ^□(□pp) 
360 ^ ^ 360 360 ^ ^ 

+ ^□2(pp) + ^□2pp + J_pn^p 
360 ^ ^ 360 360 



3360 3360 ^ 2520 ^ 

1 , ^....^ . 1 o,-l„.,_^ . 1 



1680 ' 1120 ' '^^ 1680 

+ ^^nrapp) —o'^pR + ^_pn2p 

15120 ^ ^ 3024 3780 

630^ ' ^ J 630 ^ ' ^ ^ 1260^ ^ , j 

+ —av.V^pR^" + -^V.V^PaR^" + —o(V.V^pR^n 
840 ^ 2520 840 ^ ^ ' 

+ \^—o\R^"""^R.^aB) + -^a(R^"''"^V.VaRu|3) + -^R'"""^av.VaRud 
[12600 ^ ^^^ 3150 ^ ^ ^^ 9450 m « ^^/^ 

+ -^—VaVpRnuV^VR'''^ —afVaRuuVR"'') + —^VaRuv^V^R'"^ 

12600 " ^ 6300 ^ " '^'^ > 18900 ^ 

15120 " ' 15120 75600 ^ 



50400 ' 75600 ' 37800 

^ -DR^'^nR^^ - -^—U^RR + J- □(□PP) 



75600 ^ 56700 453600 

i + 0[3ft^]. (108) 



^ -DPap + --l--n2 (pp) 



453600 129600 

Two terms in Eq. (11081) . R^'^^^'^R^^ and □P'"^nP^,^, have their numerical coefficients 
corrected from [25]. Although a2(x, x) was computed with 0[3?'^] accuracy, its expres- 
sion above is, of course, exact. 

Now we compare our results with the coefficients available in the literature. The 
main concern is the fourth coefficient a4(x, x), as expressions fll06p . fll07p for a2(x, x) 
and a3(x, x) coincide with the results obtained by other methods |lOl [291 EI]) after 
taking into account differences in definitions and curvature conventions. In [34], the 
ffist order curvature terms are presented for all Schwinger-DeWitt coefficients. This 
result is fully consistent with similar terms in coefficients fll06p - fll08l) . 

The only general result for the fourth Schwinger-DeWitt coefficient a4(x,x) be- 
longs to I. Avramidi [29]. Unfortunately, his result is not expressed in terms of an 
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irreducible basis of tensor invariants. The tedious work of the reduction to a tensor 
invariant basis, that we choose to be the basis above, was performed. It was done par- 
tially with help of tensor manipulation software packages, MathTensor [26] and Ricci 
j27j . This reduction is even more difficult than a similar procedure for JdxtTa4{x,x) 
performed in [20] because of the presence of total derivatives. For example, the re- 
sult of [29] contains a term which represents six covariant derivatives acting on the 
Ricci tensor, totally symmetrized over indeces . It was reduced to the term O^R, 
while keeping all curvature terms arising from the derivatives commuting. The use 
of computers for such calculations is unavoidable, and the ability of Ricci to perform 
index free manipulations is invaluable. The result derived this way from [29] is in a 
full agreement with (11081) . It should be noted, that we were only concerned with re- 
producing the squared curvatures terms, and the complete fourth Schwinger-DeWitt 
coefficient in a form compatible with the lower order ones remains to be derived from 
either [29] or [33] • 

The work [52j gives only pure gravitational terms, thus, only terms of (11081) in 
the square brackets can be compared with it. Our expression (11081) disagrees with 
[32] in several numerical coefficients. Taking into account, that results of [32] fail to 
reproduce correctly even the lowest order terms of the trace of the fourth Schwinger- 
DeWitt coefficient, this work should be considered erroneous. 

We also mention that it might be possible to obtain the Schwinger-Dewitt coeffi- 
cients above from the short time expansion of the trace of the heat kernel found in 
[20] by applying the variational principle f|T6l) . However, in order to derive the fourth 
coefficient a4(x,x), that expansion should have been derived up to the fifth order 
in the proper time. This fact follows directly from the Schwinger-DeWitt expansion 
ffTUDD and (dS]), 

an-i{x,x) = g~^^'^-p- / dx g^^hr a„(x,x), (109) 

and was used in the gauge field theory [23] . 

It is even easier to derive the large (sometimes called late) proper time asymptotic 
of the heat kernel. The large time expansions of the basic form factors are known 

/(-sn) = -l^ + o(l), s^oo, (110) 

Fi-sa,, -.□3) = ^(5;^ + 5^ + 5^) + O (1) , . ^ 00. (Ill) 

Their substitution to the expressions for the heat kernel's form factors leads to the 
following equation, 

k(s\x,x) = -^—g^/^i-2-P+--Ri 

-(n-)(5^)-5(4^) 
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+2- V.-P V«-7^^"^ - 2- f Vfl-7^^" V„-P 



-2I (v,i7^'^'^ v"i7^„.) + o[j?3]}. (112) 

This formula is in agreement with the result of variation of the large time asymptotic 
of the trace of the heat kernel found in 1201. 



5 Green function for the massless scalar field in 
two dimensions 

As a test of our results for the heat kernel, in this section we derive the Green function 
of the scalar field model in two dimensions. The curvature expansion for the one- 
loop effective action ([9]) in two dimensions does not exist due to infrared divergences 
appearing at every order in the curvature due to Eq. (fT5l) . The exception is the two- 
dimensional effective action for massless fields [35], which does exist and admits an 
exact form. 

How this effective action arises from the covariant perturbation theory was shown in 

Naturally, we expect this behaviour to hold also for the corresponding Green func- 
tions (fTT]) - (fT2|) . As shown in [35], we can start by deriving the exact two-dimensional 
Green function from the heat kernel in a fiat space-time (l29l) . Then we performs the 
Weyl transformation ^36j, which in two dimensions consists of only the local rescaling 
of the metric, 

^^'^(x)=e^(^)(?^'^(x), (114) 

where S(a;) is a smooth function vanishing at space-time infinity. The outcome would 
simply be, 

G{x,x) = -^g'/\x)J:ix). (115) 
47r 

This expression admits a nonlocal closed form, upon inserting the conformal factor 
S expressed in terms of the covariant metric 

^g) = ^R- (116) 

In order to arrive at this result, we proceed first by formally defining the regular- 
ized Green function via the heat kernel as, 



G''^ = / ds 







K{s) - K{s)W=o , (117) 



and then computing with the use of the explicit form of the heat kernel from the 
section [31 
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There is only one curvature in two dimensions due to the properties, 



trl = 1, 7^^, = 0, 



Therefore, the heat kernel (H6l) reduces to: 

,1/2 



oi/2 r 2 

Ki-S) = ^y- + sJ2c^9^i-Sa)R 



11 



i=l 



where the coefficients c,- and C*,- are 



(118) 



(119) 



Ci — 7:, C2 — 1, 
o 



36 



12 



Cio — 2'"'^' ~ ^'~'i'~'2' C2 — Cq — Cj — 0. 



(120) 



The final expression looks like 



K{s)-K{s)\r-_ 



9 



1/2 



4 



47r 
+ 

+ /(-sni 



l f{-sD)-l 
2 



R 



sF(-sni,-sn2,-sn3) 
□1 



2 I CI— 1 1— 1 2 |— 1 3 I ri—i 1—1 2 |— 1 3n 



+ 



2D2D 



'f{-sn,)-i\ 1 



2D2D 



(□1 - □ 



+ 



402 



:ni-n3) 



2n, 



-R1-R2 



+ 0[i?3]}. 



(121) 



A straightforward check with help of the large time asymptotics f lllOI) - flllip shows 
that the expression fll2ip behaves like ai s —>■ 00, which provides the convergence 
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of the proper time integral flll7p . To perform the proper time integration we resort 
to the technique of [91 [20], namely, we express K{s) as a total derivative over the 
proper time: 



Kis) - K{s)\r=o 



47r ds 



{n{s, □) + m^y^is, Di, Ds, Ds)) + OlR" 



where 



nis,n) = -f{-sn)R, 



(122) 
(123) 



and m'^y™ is the following function symmetrized over indices 1 and 2: 



m(s,ni,n2,n3) 



-s—^F{-sa^,-sn2,-sO, 



-□in3 + ni2_3n2n 



21-13 



11-12 



1 Di //(-sDa) /(-sDi)^ 



- Da) 



□ 1 



R1R2. (124) 



After Substituting the expansions flllOI) - (lllll) into (I122p and doing the proper time 
integration, we observe that the whole expression vanishes at large s, and the resulting 
equation is 



9 



1/2 



4:71 



niO, □) + m^y'"(0, Di, 02, Da) + 0[R% 



(125) 



As is seen from (11241) . only the linear in curvature term survives, so, the Green 
function we seek is 



An □ 



R+OIR'^ 



(126) 



Although we have only demonstrated that the second order in the curvature ex- 
pansion is exactly zero, it is clear from flll5l) - (lll6p that the answer (11261) is, in fact, 
exact, in the same fashion as the one-loop effective action (I113p is exact in two dimen- 
sions P,|2T]. We note, however, that the expression (11260 cannot be directly obtained 
from the corresponding effective action (IllSp by the generating function method [T^ 
because there is no potential term P with respect to which we can take a variational 
derivative. 

It is appropriate here to remind one how to understand and work with the results 
of this paper or any paper on the covariant perturbation theory. 



^R{x) = JdyG{x,y) R{y), 

where the Green function G{x,y) of □ is defined as, 

a^G{x,y) = 5{x,y), G{x,y)\,, 



0. 



(127) 



;i28) 
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6 Summary 



Let us summarize this study. In this work we have implemented the covariant pertur- 
bation theory for the coincidence hmit of the heat kernel. It was computed up to the 
second order in the curvatures by the generating function method from the known 
heat kernel trace, and also directly by the covariant perturbation theory algorithms. 
The basis of the second order nonlocal tensor invariants was obtained. The nonlocal 
form factors are given in the two integral representations. By taking the functional 
trace of the heat kernel, the lower, known order of the heat kernel trace is repro- 
duced. We have compared the short proper time expansion of our results with the 
Schwinger-DeWitt coefficients found in the literature. The work [29] is confirmed, 
while the other [32j is proved to be erroneuos. The diagonal Green function for the 
massless scalar field model in two dimensions was derived and shown to produce the 
closed form nonlocal expression. Due to these tests, the results of this paper can be 
deemed, via the variational principle, another verification of the third order of the 
trace of the heat kernel derived in [9l |19l [20] . 

The validity of the covariant perturbation theory obeys the condition [6], 

VV3? » (129) 

This means is the results are applicable to problems with rapidly oscillating back- 
ground fields of small magnitudes, i.e., the high frequency limit. The opposite case 
of slowly fiuctuating fields of large magnitudes, which corresponds to the effective 
potential, also presents a significant interest and was studied in ^7\ . 

The obtained results can be directly applied to the computation of the Green 
function in four dimensions. To derive it, one needs to make an additional step, 
the integration over the proper time (fTT]) . It is proved that only the first order in 
curvatures is ultraviolet divergent [101 113j, which in the dimensional regularization 
[2H] is 

G{x, xf^ = TT^-^ g'^'P, ^ - 2, (130) 

while the second and higher orders are finite and nonlocal. As was mentioned, infrared 
divergences of massless field theories do not appear at all due to the nice behavior of 
the heat kernel at large proper times (cf., Eq. f[T^ ). 

K{s)ocs-'^, s^oo, 5R^0. (131) 

which holds for any curvature order excluding the zeroth. 

The found nonlocal contributions of the Green functions are responsible for a 
variety of physical effects [18]. All physical information is contained in the nonlocal 
form factors. However, one needs to transform these form factors to the massive 
Green functions [3 [6l [TJ [8] . The first order form factors turn, upon the proper time 
integration, to 

-'■'(^) = r*"'(;^-;;?T7?)' ''''' 
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where /x^ is the parameter of uhraviolet regularization. For the second order form 
factors, the following construction can be derived [7], 



These spectral forms encode the information about the filed model into a set of 
spectral densities p{mi^ , nii^ , rris^) [39]. Another step towards physical applications 
is going from the Euclidean signature metric to the Minkowskian one. The unique 
Green function of the Euclidean space gets replaced by one of many. Which one is 
being defined by the physical problem under consideration via the integration contours 
|lUj . If one is solving a Cauchy problem with data in the remote past, for expectation 
values of fields in the standard m-vacuum, this procedure boils down to using the 
retarded massive Green functions in the effective equations [5l H8j. 

The very first application of the nonlocal effective action was quantum electrody- 
namics where the correct derivation of the electron magnetic moment, free of infrared 
divergences, was done [H]. General results of the covariant perturbation theory are 
quite universal and applicable to a wide class of physically interesting quantum field 
theories that are specified by the differential operator -F(V) of the form ([3]). Starting 
from the model's Lagrangian one can write down -F(V) and explicit forms for the 
curvatures P and TZ^^. The prescriptions on reducing various field models to the 
standard form can be found in [13]. 

For the example of quantum electrodynamics coupled to gravity, the whole struc- 
ture of the nonlocal effective action was already worked out by G. Shore [12]. It 
seems instructive to show here explicit forms for the CPT curvatures of this model, 
omitting the model's description and explanations that can be found in [12]. The po- 
tential term is P = iea^'^F^y — ^^-Rl, and the commutator curvature takes the form 
of it^y = ieF^yl + \cr^^R^u\p) where standard QED notations for the electromagnetic 
tensor and the spinor matrices are used. We should emphasize that although CPT 
was developed for massless field theories, its generalization for massive theories like 
QED is simple. It boils down to an additional factor e"*'" under the proper time 
integral [13] that involves only form factors but not the curvatures. Explicit examples 
of such calculations can be seen in [121 113] . 

To look at the relation of the used techniques and their results to other approaches 
in quantum field theory, we mention the so-called string inspired or world line methods 
that were developed in recent years, e.g., see review |46j. Their core is the world 
line method [15]. It is a path integral perturbation theory where the world line 
parameter is effectively the Schwinger's proper time [16]. Indeed, as was shown in 
j47j the form factors of the one-loop effective action, obtained by the world line 
method coincide with those calculated in the covariant perturbation theory ^ [21] , 
in flat spacetime. Nevertheless, these approaches are different, and each has its own 
domain of applicability and the set of problem it is best suited for. 

We expect that our results will be applied to gauge field theories. The chiral 
perturbation theory was one of the applications of the Schwinger-DeWitt series [23] . 
so, more powerful techniques must find their use as well. Still, the main motivation for 



/ 



(imi^(im2^(im3^p(mi, m2, ^3) 



(133) 



'0 (mi2 - □)(m22 - □)(m32 - □)■ 
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the development of the covariant perturbation theory was quantum gravity problems 
[H [IHl [39] . Pursuing the systematic development of the effective action method for 
particle creation [U] G.A. Vilkovisky has solved the quantum gravitational collapse 
problem that appeared to have an intriguing and unexpected solution [49J. 

Further progress in the development of the covariant perturbation theory may lie 
in the construction of a partial summation of the curvature expansion and/or the 
derivation of the heat kernel with separated points, K{s\x, y), but such developments 
are beyond the scope of this paper. 



Acknowledgments 

We acknowledge the support by NSERC of Canada at early stages of this work. 
The author is very grateful to A.O. Barvinsky and G.A. Vilkovisky for the long 
collaboration on the the covariant perturbation theory. 



Appendix A: The reduction equations for the form 
factors 

This section reproduces the equations for the a -polynomial reduction of the form 
factors that were derived in [HI [9] . 

By eliminating all polynomials in a -parameters in the form factors of the heat 
kernel, the form factors can be explicitly expressed through the basic form factors. 



= f^dae-"^^-"^^, (134) 
Jo 

i^(6,6,e3) = (e^^)3 = 

/ (fa 5(1 - ai - ^2 - as) exp(-a;iQ;2'C3 - ^20(3^1 - 0103^2)- (135) 

Ja>0 

This kind of form factor representation is called, therefore, the explicit representation. 

After the use of the delta-function in fll34p and fll35p . there are two types of 
Q!-monomials: 



daa^exp — a;(l — a)^ 



(136) 



For the case 



n „ m^sQ 
tl«2 



1 /■!— 02 

da2 / dai a\a\ 
h 



tti e 

X exp — «2(1 — «! — a2)^i — 0L\{\ — ai — 02)^2 — Oiia2$,3 
integration by parts produces the following equation: 





1 r 

da —a"^ exp — a{l — a)^ 
eta L 



n = 
n > 



(137) 



(138) 
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which yields the recurrence relations, 



«ie-°^°^«)^ = \ (e-"^"^«)^ , (139) 



-2(^-l)(«r'(^ ^ j^^, ^>2. (140) 

They allow one to express all integrals (11361) through the basic form factor (11341) . Since 
the recurrence relations imply a division by the appearing subtractions maintain the 
analyticity of the integral (I136P in ^ at ^ = 0. For the form factor with subtractions 



in terms of (11341) one has 

lp:^\^^^m_i^ ,141) 

Elimination of the polynomials in a from the form factors (11371) is based on inte- 
gration by parts: 

rl |'l-«2 d , 

\ ^ V y/2' Q42) 

(<a^e-"i"2«3^2^ ^>0. 

A similar expression involving the differentiation with respect to a2 is obtained from 
(11421) by transmutations of the indices 1-2. The first order form factors appearing on 
the right hand sides of them are subject to the recurrence relations above. By perform- 
ing the differentiations on the left hand sides of these two formulas, one obtains two 
linear algebraic equations for the quantities (^ai~^^a^e^^^^ , (^a"a™'''^e*^^^containing 
the highest order monomials. The discriminant of this linear system is 

A = + e2 + el - 2^16 - 26^3 - 2^6, (143) 

and the recurrence relations obtained in this way are of the form 

6(^3 + 6-6) / 



I ^ (1^3 - 6 - 6) / „ m-l sQ 

+ m (a, an e , 

A \ ^ /3 

(6 + ^1-6 



A 



(<a^e-'^^"^«^)^ + /3(n,m), (144) 
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P{n, m) = 0, n > 0, m > 



(145) 



p{n, 0) = - 6 - ^i) ^^"^--1-26^^ ^ ^ > (146) 
/3(0,m) =2^(a7^e-"i"2«i)^, m>0 (147) 

/?(0, 0) = 2^ (Je-"^"2«i)^ + - 6 - 6) ^g-aia26^^ _ (148) 

An expression for <^a"a^'''^e*^^^ can be obtained from fll44l) - fll48p by the 1-2 indexes 
transmutation. Together with fll39l) - fll40p these relations make it possible to express 
all integrals (11371) through the basic form factors /(^) and F(^i, ,^2, ^s)- For the 
combination fl65l) with a subtraction, one has 

e^^- 1)^ = ^(^1,6,^3)-^, (149) 

in terms of (11351) . Analyticity in 6, 6, ^ holds now only in the sum of the form factors 
on the right hand side of (I144p . and it is a nontrivial fact that, when these form factors 
are expanded in power series in ^ the denominator A gets always cancelled. The 
mechanism of maintaining analyticity is based on the existence of linear differential 
equations which the functions (11341) and (I135p satisfy. 

The differential equations for the basic form factors can be derived with the aid of 
the recurrence relations above. From (11341) and (11400 one has the following equation 
for the function f{^): 

-> - \m.\^. (150) 

Similarly, one obtains the equation for the form factor (I135p : 

-^i^(6,6,e3) = ^[(6 -6 -6) A 

+ 6^3(26^3 - ^2 - + e?)] 6, ^3) 

1 8666 + (6 + 6-^1) A ,,, , 
+ 2 /(^i) 

, p,66(6 - 6 - 6) r-f(- N 
+ 2 / j 

^2 ^26(6-6-6) ^(^^)_ (151) 

Finally, as a consequence of these equations, one can derive an equation for the form 
factor (|T35l) function of s: 

-s— F(-sni,-sn2,-sn3) = [s ' ^ ^ + ij f(-sDi, -sD2, -503) 
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2D jy 'I 

+ ^/(-^Da), (152) 

D = Di^ + + _ 2nin2 - 2nin3 - 2n2n3. (153) 
This is an equation used for the study of a two-dimensional field model in section [51 

Appendix B: The explicit representation for the sec- 
ond order form factors 

The second order form factors of the heat kernel are expressed through basic form 
factors (HMD, (USSI) and A ffHH]) : 



G'i(ei,6,e3) = i^(ei,6,6 



(154) 



G'2(ei,6,e3) = i^(ei,6,6 



A2 



(6 + 6-6)(6 + 6-6) 



+ 1(6 - 6 - 6)] - /(6)^(6 + 6 - 6) 
-/(6)^(6 + 6-6) 
-/(6)^(-6666-366' 

- 366' + 366' + 366' + 6' + 66' + 66' - 6' - 6'), 



(155) 



G3(6,6,6) = i^(6,6,6 



+ 46^66 + 6' + 6' + 6' - 66'6' + 266' 



3A2 

2f 2 



(-2666' - 2666' 



266 
4^ 



66'6 



+ 266' + 266' + 66'6' - 46'6 - 466') - ^ 

+ /(6)f^(266 - 6' + 66 + 66 - 6') 

+ /(6)^(6 + 6 - 6)(6' + 266 - 6' - 6') 

+ /(6)^(2666 - 66' + 66' + 6' - 66' 
+ 6' - 66' + 366' - 366' - 6'), 



(156) 



G4(6,6,6) = i^(6,6,6 



:(6'-466' 



- 166'6' 



36A4~ 

686'6' - 1006'6' + 686'6' 



166'6' 



26 



- 46^6 + 26' + 16666' - 6066'6' + 866'6' 

+ 6866'6' - 4866'6' - 4666' + 866' + 966'6'6' 

- 1366'6'6' + 366'6'6 - 166'6' + 646'6'6' - 2866'6' 
-406'6' + 466V) 



- 66' - 6' + a' + 1066'6 + 5^2%' - 6' - ^^2%' 
+ 566' + 566' + 236'6'6 + 5766'6' + 126'6'6 

- 296'6'6 - 206^66' - 6666' + 3766'6' 

+ 116'66' + 366' - 426'6'6' - 6866^6' - 2266^6 
+ 146'6'6' - 146'6'6' + 366'6' - 56'6' + 56'6' 

- 366' + 56'6' - 96'6' - 276'6' + 276^6' 

- 56'6' + 56'6') 

- /(6)^(-6' - 266' + 386^6' - 906^6' 

+ 906'6' - 386'6' + 266' + 26' - 1066'6 + 506'6'6 

- 246''6'6 + 26'^/6 - 2066^6 + 146'6 - 66626'6' 
+ 366'6'6' + 6266'6' - 546'6' - 446'66' + 386'6') 

+ {^^^) §(46^6 + 36^ + 26^6^ - 866^ 

- 6' - 866' + 466' + 16666' - 66'6' + 66'6' 
-6' -8666' + 466') 

^ (^^^^6^) ii^^^^'' " ~ ^°^^'^'' 

+ 526'6' - 2666' + 26' + 3666'6 - 526'6'6 

+ 866'6 - 66'6 + 186'66' + 46'6') , (157) 

^5(6, 6, 6) = (i^(6,6,6) - ^) ^ - /(6) ^~g|^^^'^ 



+ ^(36' -1666' + 46 V 

+ 246'6' - 2666' + 86' + 2866'6 - 526'6'6 
+ 126'6 + 266'66' - 206'6')^ 

+ (i^(6, 6, 6) - ^) ^(6' - 466 + 26' + 466) 




(158) 



G6(6,6,6) = i^(6,6,6) 



(-6 + 6-6) X 



27 



X (-6 - 6 + e3)(-6 + 6 + 6) + ^ 



+ + 6 - - 6 + 6) 

- - 6 + e3)(-a + 6 + ^3) 

-/(6)^(ei'-2ei6 + 6'-C3'), 



(159) 



^7(6, 6, 6) = 6, e3)|(-e3 - 6 + 6) - /(6)^(-e3 + 6 + Ci) 
+ /(6)|-/(e3)^(e3 + 6-ei) 



(160) 



46 



^8(6, 6, ^3) = i^(ei, 6, ^3) - ^(-^3 - 6 + 6)' - ^ 



+ /(6) 



[-Cs + 6 + ei)(-46e3 - 466 + 36' - 266 + 6' + 6') 



86 



+ /(6)^(6 + 6-6) 
+ /(6)^(6 + 6 - 6)(-466 - 466 + 36' - 266 + 6' + 6') 

-(e7^(/(^0-/(6))|:, 



6 



(161) 



G9(6, 6, 6) = i^(6, 6, 6)^(-26' + 266 + 6') 

-(F(6,6,6)-^)^(26-6) 
/(6) - 1 326 



/(6) 



266 



+ 



6 + 6-6) 



6 

(26^-86'6 + 66V-i66'6 



/(6) - 1 1 /o,- 4 ot St , fit 2t 2 ict3 



6 A2ei6 

+ 166'66 + 366'6' - 20666' - 3266' + 56'), 



(162) 



:(6 + 6-6)'(-2666' 



^10(6, 6, 6) — -^(6, 6, 6) 2^ ^4 ' 

- 2666' + 466'6 + 6' + 6' + 6' - 66'6' + 266'' + 266' 

- 66'6' + 266' + 266' + 66'6' - 46'6 - 466') 

+ ^(-14666' + 1466' + 4066'6 - 14666' 

- 2266' - 126'6' - 46'6 - 466' + 66'6' + 6' 
+ 196' - 126'6' + 1466' - 2266' + 6')" 



28 



+ (^^(6,6,6)-^) 



486 

+ /(6)7^(6 + 6 - 6)(-i666'6' - 2066'6' + 24666' 

+ 866' - 646'6' + 356'6' + 376'6' + 866' - 646'6' 
+ 376'6' + 356'6' - 866' - 866' - 96' + 6' + 6' 

- 206^66' - 346'6'6' - 6666' - 1666'6' + 24666' 

- 66'6 - 666' - 206'6' + 156'6' + 156'6') 

- /(6)^(6 + 6 - 6)(-2666' - 2666' + 466^6 

+ 6' + 6' + 6' - 66'6' + 266' + 266' - 66'6' 
+ 266' + 266' + 66'6' - 46'6 - 466') 

- /(6)7Ai(-4i66'6' + 346'6'6' + 306'6'6' + 416^66' 

66 A-* 

- 4466'6' + 2066'6' - 26666' + 96' + 6' + 8466'6' 

- 8166'6' + 38666' + 766' + 336^66' + 2566'6' 

+ 26'6'6' - 5566'6' - 436'6' + 276'6' + 296'6' - 296'6' 

- 6' + 6666' + 66' - 766' + 216'6' + 766' 

+ 356'6' - 356'6' - 2i6'6' - 1766' + 996'6' - ioi6'6' 

- 276'6' + 456^6' - 966') 

+ [^^^) ^(^2 + 6 - 6)(-2666' + 10666' 

+ 40^i6'6 + 36' + ^i' + 6' + 126'6' - 1066' - 666' 
+ 446^6' - 266' - 4666' + 66'6' - 46'6 - 466') 

- [^^^) ^(-466 + 26' + 66 + 66 + 26' - 36') 
_ ^ ms)- 1 ^ _L_(4^2^^^33 + 18^36^6' + 20666' - 36' 

- 226^6' + 186'6' + 1366' - 766' + 6' - 426'6' 

- 76666' + 666'6' + 4366' + 26'6' - 1866'6' + 66' 

- 566' - 6' - 106'6' + 106'6' + 566') 

-7^^(/(6)-/(6))7^ 

(6 - 6) ^ ^ 66 

_ 1 / /(6)-i _ /(6)-i \ 1 
(6-6) V 6 6 7 6' 

Gil (6, 6, 6) = i^(6, 6, 6) [^(26' - 86'6 + 66'6' - 46'6' 
1 fi 

+ 4666' + 6') - ^(-36' + 36'6 + 46'6 



29 



+ (i^(6,6,6^ 



■(26' + 466 -466 + 6 



- /(ei)^%^(-6 + 6 - 6)'(-6 + 6 + 6) 



/(6 



A4 

1 



-(-26' + 186'6 - 426^6' + 266'6' 



2A466 

+ 146'6 - 766^66 + ir86'6'6 - ii66'6'6 
- 426^6' + ii06'66' - 686'6'6' + 706'6' 



406'66' - 306'6'6' 



706'6' - 506^66' 



+ 426V + 26666' 



1466' + 6 



+ 



7(6) -i \ 326 
7(6) -i\ 1 



(36' - 66 - 26' - 66 + 466 - 26 



-26' + 66'6 - 46'6' + 186'6 



6 J A366 
'66 - 426'6^ 
+ 686'6' - 20666' - 4266' + 56') 



+ 246^66 - 426'6'6 - 526^6' + 526'66' 



(164) 
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